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O ' Abstract 

o. 

Size-selective concentration of particles in a weakly turbulent protoplanetary nebula 
may be responsible for the initial collection of chondrules and other constituents into 
Q-r primitive body precursors. This paper presents the main elements of this process of 

ry\ . turbulent concentration. In the terrestrial planet region, both the characteristic size and 

size distribution of chondrules are explained. "Fluffier" particles would be concentrated 
in nebula regions which were at a lower gas density and/or more intensely turbulent. 
The spatial distribution of concentrated particle density obeys multifractal scaling, 
suggesting a close tie to the turbulent cascade process. This scaling behavior allows 
predictions of the probability distributions for concentration in the protoplanetary 
nebula to be made. Large concentration factors (> 10 5 ) are readily obtained, implying 
that numerous zones of particle density significantly exceeding the gas density could 
exist. If most of the available solids were actually in chondrule sized particles, the 
ensuing particle mass density would become so large that the feedback effects on gas 
turbulence due to mass loading could no longer be neglected. This paper describes 
the process, presenting its basic elements and some implications, without including the 
effects of mass loading. 
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1 Background and Introduction: 



>< 



Primitive (unmelted) chondritic meteorites are composed in large part of mm-sized, once- 



molten silicate particles (chondrules) and metallic grains out of mineralogical equilibrium 
with each other. Many chondrites contain inclusions of refractory minerals that have been 
dated as the oldest objects formed in the solar system (MacPherson et al 1989, 1995). The 
chemical, isotopic, mineralogical and petrographic properties of individual chondrules them- 
selves imply that independent entities were melted by some "flash heating" event in the 
gaseous protoplanetary nebula, and remained molten for fairly short times (less than an 
hour; Jones et al 2000). Chondrules are diverse in chemistry but are narrowly size-sorted 
(Grossman et al 1989, Jones and Brearley 1999), apparently by their aerodynamic cross sec- 
tion (Dodd 1976, Skinner and Leenhouts 1991, Keubler et al 1999); the least mechanically 
evolved pieces in chondrites have the appearance of being gently brought together (Metzler 
et al 1992) and subsequently compacted to solid density. Various hypotheses have been ad- 
vanced to explain these properties (see, e.g., Boss 1996, Hewins et al. 1996, Hewins 1997, 
Connolly and Love 1998, and Jones et al 2000 for recent general reviews and discussion), 
but the puzzle remains unsolved. The fact that many primitive meteorites are composed 
of up to 70-80% chondrules by volume (Grossman et al 1989) implies that the chondrule 



formation and accumulation processes were of widespread and significant importance in the 
very earliest stages of the accretion of the asteroidal objects which provide the parent bodies 
for these primitive meteorites. Since the terrestrial planets apparently formed from various 
proportions of known meteorite types, one expects that the process extended beyond the cur- 
rent asteroid belt and that understanding these early stages is important for understanding 
planetary accretion overall. These stages of circumstellar disk evolution relates to regions 
and epochs when the particulate opacity is high, but when some particle growth has occurred 
- and are thus also of great interest for study at infrared and millimeter wavelengths. 

Prior work (Dubrulle et al 1995, Cuzzi et al 1996) has pointed out that, unless the tur- 
bulent kinetic energy in the nebula gas is vanishingly small (see Section 2), chondrule-sized 
particles are unable to settle individually to the nebula midplane, where most growth to 
planetesimal size must occur. Instead, Cuzzi et al (1996; henceforth CDH96) proposed that, 
following their initial melting, and throughout multiple recurrences of similar heating events 
(Wasson 1996, Connolly and Love 1998, Desch and Cuzzi 2000), chondrules pursue an ex- 
tended free-floating existence under plausible conditions of nebula gas density and turbulent 
intensity in the terrestrial planet region, successively encountering zones of varying concen- 
tration enhancement (section 5.1), until by chance they encounter an unusually dense zone 
where they might physically coalesce into much more massive, but still not solid, entities. In 
a second stage, such entities - dense clusters of particles - might have enough coherence to 
resist disruption as they settle to the midplane, or to be collected in a different process into 
the cores of the largest eddies (discussed below) for subsequent accumulation into planetesi- 
mals. Or, dense zones may only provide environments of enhanced collisional accumulation 
of chondrules. The second stage remains unstudied and qualitative; here we focus on the 
first stage. 

The process of interest is known as preferential concentration or turbulent concentration 
(TC). The nature of this process is that isotropic, homogeneous, 3D turbulence provides 
numerous fluid zones of low vorticity and high strain in which particles having a narrowly 
defined range of aerodynamic properties can be significantly, but usually briefly, concentrated 
(discussed further below). This somewhat counterintuitive effect was first alluded to the- 
oretically (Maxey 1987), subsequently demonstrated numerically (Squires and Eaton 1990, 
1991; Wang and Maxey 1993), and recently demonstrated experimentally as well (Fessler 
et al 1994; for a review see Eaton and Fessler 1994). The concentration factor C is the 
ratio of the local particle density to its global average. In numerical studies to date (Squires 
and Eaton 1990, 1991; CDH96, Hogan et al 1999), turbulence with Reynolds numbers of 
10 2 — 10 3 contains dense zones where the concentration factors C reach 40-300 for optimally 
concentrated particles. We extended numerical studies of this effect, and applied scaling 
relationships to predict its behavior under protoplanetary nebula conditions (CDH96, Cuzzi 
et al 1998). We found that under canonical inner nebula conditions, particles with the size 
and density of chondrules would be optimally concentrated. Here, we show that the process 
is easily generalized to a wide range of fluffier particles in lower density regions of the nebula 
{section 3). 

We note that TC is a quite different process than the superficially similar effect in which 
far larger (meter-and-larger-radius) particles with stopping times comparable to or larger 
than the eddy times of the largest eddies (with frequencies comparable to the orbit period) 
can accumulate near the centers of such eddies (Barge and Sommeria 1995, Tanga et al 1996, 



Bracco et al 1999). Klahr and Henning (1997) showed concentration of mm-sized particles 
within large, slow, 2D circulation patterns (not turbulent eddies). However, unless these slow 
circulation patterns represent the primary kinetic energy reservoir of the (non-Keplerian) 
fluid motions, concentration of such small particles will not occur. That is, if the nebula 
turbulence has a normal 3D cascade with its energy peak at spatial and temporal scales which 
are a small fraction of the nebula scale height and orbital period (as normally implied by 
turbulent "a-models" described in section 2), chondrule-sized particles will diffuse faster than 
they can be concentrated by such large, slow circulation patterns. For example, numerical 
calculations by Supulver and Lin (1999), which modeled eddy motions on a wide "inertial 
range" of spatial and temporal scales, as well as our own calculations, show that chondrule- 
sized particles do not accumulate in large eddies for two reasons: first, their trajectories are 
mixed by the smaller eddies, and second, "real" eddies in homogeneous turbulence - even the 
largest ones - dissipate within a single overturn time. Under these circumstances (assumed 
in this paper), mm-sized particles would be nearly uniformly dispersed by 3D turbulence, as 
one would naively expect for such small particles which are trapped to nearly all fluid parcels, 
if it were not for the role of turbulent concentration as described herein. The true energy 
spectrum and dimensionality of nebula turbulence is, however, not currently understood, 
and is the subject of several active research tasks. 

Our emphasis in testing these concepts has been to compare our predictions with the 
properties of easily identified and studied chondrules and chondrites. For instance, we have 
now studied the size distribution of preferentially concentrated particles in detail, and found 
it to be insensitive to, or independent of, Reynolds number (Hogan and Cuzzi 2000); here we 
show that this predicted size distribution is in very good agreement with a typical chondrule 
size distribution (section 4)- 

While C does increase systematically with increasing Reynolds number (CDH96), our 
prior estimates of concentrations at nebula Reynolds numbers, which are plausibly far larger 
than those accessible to numerical modeling, had required sizeable extrapolations. We have 
more recently shown that the spatial structure of the concentrated particle density field is 
a multifractal which has Reynolds-number-independent properties (Hogan et al 1999) and 
here will use this result to provide a firmer basis for predictions under nebula conditions 
(section 5). 

The evolution of extremely dense clumps, within which interparticle collisons might en- 
trap particles (CDH96) remains unstudied, and will require a better understanding of the 
behavior of turbulent concentration when the particle mass density exceeds that of the gas, 
and of particle ensembles whose density is large enough to affect the gas flow properties. In 
section 6 we discuss this important effect. 

2 Turbulence and turbulent concentration (TC): 

Homogeneous, isotropic, 3D turbulence is characterized by a cascade of energy through a 
range of scales, known as the inertial range, from the largest (or integral) spatial scale L, 
having associated velocity Vl, to the smallest (or Kolmogorov) scale r) where it is dissipated 
(Tennekes and Lumley 1972, Hinze 1975). The intensity of the turbulence is characterized 
by the Reynolds number, which can be written Re = (L/r/) 4 / 3 . Since dissipation occurs 



primarily on the small scales where molecular viscosity comes into play, more energetic 
(higher Re) flows can drive turbulence through a wider inertial range, to smaller 77, for any 
given viscosity. 

Re is ordinarily defined as Re = LV^jv m ^ where v m is the molecular viscosity; Re here is 
the ratio of transport by macroscopic motions to that by molecular motions. This definition 
combines the velocity and length scales into a turbulent viscosity vt- The Reynolds stresses 
which appear in angular momentum transport equations are often modeled by this sort 
of turbulent viscosity (the "alpha model" of Shakura and Sunyaev 1973). However, in a 
Keplerian disk, the mere existence of turbulent motions (turbulent kinetic energy, leading 
to diffusivity of scalars) does not necessarily imply viscous transport of angular momentum 
(by a Reynolds stress, or torque, acting as a positive turbulent viscosity; Prinn 1990). This 
distinction is related to the possibility that the on-diagonal and off-diagonal terms of the 
stress tensor might have very different relative strengths in Keplerian disks than in more 
familiar turbulent environments (Kato and Yoshizawa 1997). While astrophysical "alpha- 
models" of protoplanetary disks emphasize turbulent viscosity Vt in its angular momentum 
transport role, TC is more closely related to scalar diffusivity, or turbulent kinetic energy 
per unit mass density k. Therefore, we distinguish two types of the familiar astrophysical 
a. We first discuss the familiar Shakura- Sunyaev prescription, which defines the turbulent 
viscosity as vt = a v cH , where c is the sound speed and H is the vertical scale height. 

Without external drivers, the overall nebula is likely to be in a regime of Rossby number 
Ro = Ql/Qq of order unity, where VLl is the largest eddy frequency and Qq is the orbital 
frequency. We thus rewrite vt in terms of these fundamental properties: 

Vt = LV L = L 2 n L (=Vl/n L ) = ai/ cH = a u H 2 n (=a u c 2 /n ); (1) 

Thus 

/n„V' 2 „ _/nt\ 1/2 



L = H ^\t) • v ^ c ^\JTj ■ (2) 

With the stipulation that the largest truly turbulent eddies (those that participate in the 
turbulent cascade) have frequencies Ql no smaller than, but probably comparable to, the 
local orbital frequency Qq, Vl = c^fcT^ and L = H^fcT^. Occasionally it is assumed that 
L rs H, but this implies Vl = v T /L = ca u and thus VLl = Vl/ L = ca v /H = a v Q << Q , an 
implausible situation for genuine turbulence. The distinction is important for us, as we need 
to scale turbulent velocities and lengthscales independently. The important point is that 
the scaling parameter a v is "shared" by the length and velocity scales, rather than being 
associated with one or the other. 

A scaling analysis cannot establish whether fluid motions with arbitrary L and Vl do 
in fact provide a positive Reynolds stress or "turbulent viscosity" given by their product; 
conversely, estimates of the magnitude of Vl and L from an observed turbulent viscosity 
vt an d associated a u , as above, may misrepresent their actual magnitudes. This is because 
certain types of spatial correlations between "random" fluid motions are needed to provide 
a positive Reynolds stress, or turbulent viscosity, and this may not occur in systems which 
are strongly influenced by rotation or otherwise strongly perturbed (Prinn 1990, Kato and 
Yoshizawa 1997). 

Another approach to determining Vl is based on the turbulent kinetic energy per unit 
mass k and an associated a^: k = |V^ = \cikC 2 , or Vl = cJat/S. The mere presence 



of turbulent fluid motions with k = akC 2 /2 = 3V£/2 is sufficient to produce turbulent 
concentration, with no additional uncertainties about the degree and sign of the correlation 
between orthogonal components of the fluid motions as in a v (Prinn 1990). While this 
definition provides no insight into the turbulent length scale L, we presume by analogy to 
the above argument that L = H^fa^. With regard to the nebula, where properties are 
uncertain, we suppress the factor of v3 (effectively ignoring the distinction between a single 
component of Vl and its magnitude) and approximate the turbulent Reynolds number of 
the nebula by Re = a^cHjv m . Angular momentum may still be transported by a turbulent 
viscosity vt = a v cH , but whether a v ~ at or not is peripheral to this work (Kato and 
Yoshizawa (1997) find that a u < «&). 

Possible sources of k include forcing by ongoing infall onto the disk in the very early stages 
(Cameron 1978, Prinn 1990), turbulent convection powered by release of gravitational energy 
into heat as the disk evolves (Lin and Papaloizou 1985, Cabot et al 1987; Goldman and 
Wandel 1994; Bell et al 1997), enforced Keplerian differential rotation (Dubrulle 1993), and 
magnetorotational instability or MR1 (Balbus and Hawley 1991, 1998). Typical estimates of 
ojfc from these sources are ~ 10~ 4 - 10 _1 , and thus Re = 10 8 - 10 11 . 

However, uncertainties remain with all of the above mechanisms. The infall stage lasts 
only a relatively short fraction of the evolution lifetime of typical nebulae, and the magne- 
torotational instability will not act in the dense inner scale height of the disk where most of 
the mass and chondrules probably reside (Gammie 1996) and perhaps not even high in the 
nebula (Desch 2000). The validity of convection and differential rotation, acting by them- 
selves, rests on their uncertain ability to be self-sustaining. This requires turbulence to be 
able to transport sufficient angular momentum outwards that mass can evolve inwards, re- 
leasing gravitational energy to be converted into turbulence. Several studies (most recently 
Stone and Balbus 1996) found that convective turbulence fails to produce outward angu- 
lar momentum transport (positive Reynolds stresses). However, some recent 3D numerical 
studies, which better capture large azimuthal structures, seem to imply that it may in fact 
be able to do so (Klahr 2000a,b). Of course, as grains accumulate, opacity decreases and 
thermal convection may weaken. Differential rotation is free of this limitation (Dubrulle 
1993), but numerical and analytical arguments by Balbus et al (1996) question differential 
rotation as a source for turbulence based on the energetics and stability of Keplerian disks. 
However, Richard and Zahn (1999) have suggested, based on laboratory analogues, that in- 
stability to turbulence in such systems requires a higher Reynolds number than accessible to 
current numerical models. Kato and Yoshizawa (1997) have shown that Keplerian rotation 
probably does not preclude some true turbulent viscosity (ie. positive Reynolds stress, or 
a v > 0), but as a much smaller fraction of the ambient turbulent kinetic energy (at) than 
under non-Keplerian conditions. Thus, it also remains in doubt whether turbulence driven 
by differential rotation alone can be self-sustaining. Non-steady situations might even need 
to be considered. 

In short, the source of angular momentum transport is poorly understood, and it is not 
clear how disks evolve at all. However, while we still don't understand the mechanism that 
allows them to do so, protoplanetary disks are observed to evolve with mass accretion rates 
in the range of M = 10~ 7 — 10~ 9 M Q /yr (Hartmann et al 1998) and thus have an associated 
gravitational energy release rate per unit area of E grav = 3GM*M/4nR 3 (Lynden-Bell and 



Pringle 1974). Here, G is the gravitational constant, M* is the mass of the central star, and 
R is the distance from the central star. 

Turbulent kinetic energy is dissipated at a rate of approximately E k = 2kQoPgH = kQoIl, 
where p g is the local gas density and £ = 2p g H. If the gravitational energy E gTSuV is released 
where most of the mass resides, and converted into mechanical turbulence with efficiency £r, 
then for a disk in steady state with M independent of radius, -Eg rav ^T — Ej,, or 

AnR 3 = ^oS = —j—. (3) 

Substituting fi 2 , = GM*/R 3 and c = #fi , 

3fe ... 

" fc = 2^S- (4) 

For a typical "canonical" nebula with H(R) = 0.05AU(i?/lAU) 5 / 4 , c(R) = LOxlO^iZ/lAU)- 1 / 4 
cm sec -1 , and S(i?) = 1700jF(i?/lAU)~ 3 / 2 g cm -2 , with T being any enhancement in mass 
over the "minimum mass" nebula density (Hayashi 1981; see also Cuzzi et al 1993), then 

o / R \ 1/2 M ?t 

°' sUxr U) Wm-^T- (5) 

Some recent numerical calculations show £t to be at least several percent (H. Klahr, personal 
communication 1999); thus a k ~ 10~ 3 — 10~ 4 seems not to be out of the question in the 
2-3 AU region where meteorite parent bodies form. Conversely, Dubrulle et al (1995) have 
shown that, in order for chondrule-sized particles to settle into a midplane layer having a 
density approaching that of the gas (a layer of thickness 10~ 2 — 10~ 3 H), a^ would need to 
be in the 10~ 8 — 10~ 10 range (see also CDH96). For a k to be this low, equations 4 and 5 
show that the conversion efficiency into turbulence would need to be less than 10 -5 for a 
disk with M = lO _8 M /yr. Here we simply presume the presence of nebula turbulence at a 
weak level and explore the consequences. Henceforth, we identify the nebula "a" as a k , and 
treat Re and its associated a as characterized by L and Vj,. 

We assume the turbulence has a Kolmogorov-type inertial range, within which each length 
scale / is characterized by velocity v(l) = Vl^I/L) 1 ^ 3 (Tennekes and Lumley 1972, CDH96) 
and eddy frequency u(l) = v(l)/l = Q (l/L)~ 2 / 3 , where the frequency of the largest eddy 
oo(L) = Vl/L is set equal to Q , the local orbital frequency. Because the most interesting 
scales for particle concentrations are on the order of rj = LRe~ 3 / A « L, and even L ps 
H^fa « H, deviations from isotropy due to rotation are not a major concern. Particles 
smaller than the gas mean free path (i.e., smaller than several cm radius under nebula 
conditions) are in the Epstein drag regime (Weidenschilling 1977), and have a stopping time 
due to gas drag which is 

t s = rp s /cp g , (6) 

where r and p s are particle radius and internal density. The particle Stokes number Sti = 
t s u(l) determines the particle response to eddies of a particular size and frequency; previous 
studies have shown that the optimally concentrated particles have St v = t s u(r]) ~ 1 (Eaton 
and Fessler 1994; Hogan and Cuzzi 2000), that is, their stopping time is comparable to the 
Kolmogorov eddy turnover time. 

6 



3 Generality of turbulent concentration 

Solving the relation St v = t s u(rj) = 1 under nebula conditions, CDH96 concluded that the 
optimally concentrated particles in the terrestrial planet region of a minimum mass nebula 
with a ~ 10 -4 — 10~ 3 would have radius and density comparable to those of chondrules. 
Here we generalize TC to the entire range of nebula conditions. Using the definition for t s 
(equation 6) to rewrite the expression for St v = t s uj{rj) = 1, 

rp = t =cpt = c ii = c p 9 = c pyi 2 m 

CDH96 Substituted v m = cA/2 = m H2 c/2p g <r H2 , where A is the molecular mean free path, 
and m,H 2 = 3.2 x 10~ 24 g and a H2 = 5.7 x 10~ 16 cm 2 are the mass and cross section of a 
hydrogen molecule. However, their expression for A did not include the finite size or the 
Maxwellian velocity distribution of the gas molecules (cf. Kennard 1938; equations 106d 
and 126b), so their A was too large by a factor of 4V2. Rearranging terms from equation 
(7) and correcting this oversight, we find 

-fe) , "(r-"^(f) v, (^r«-. <*> 

where S(i2) = 2p g H is surface mass density at some distance in the nebula. In the last 
term we have adopted a canonical radial dependence of S(i?) = 1700jF(i?/lAU) _3//2 for a 
"minimum mass" nebula, with JF being some mass enhancement factor. 

This relationship is shown in figure 1 for several different typical locations and two 
different nebula masses. The range of rp s for most chondrules (Ch), obtained using data in 
Grossman et al (1989), is mapped along the line for R = 2.5 AU, and indicates the range 
of a required to concentrate them selectively. Higher values of a (more intense turbulence) 
select smaller and/or lower density particles. In both cases, porous aggregates (PA), having 
considerably lower radius-density product than chondrules, can be optimally concentrated 
at the low gas densities which characterize the outer planet region or low density regions 
high above the nebula midplane. Such porous objects are easily produced because the low 
relative velocities in turbulence of both porous, low-density aggregates, and their constituent 
monomers, lead to large sticking efficiency and minimal disruption (Weidenschilling and 
Cuzzi 1993; Chokshi et al 1993; Dominik and Tielens 1997). It has also been suggested 
that extremely low density, intensely turbulent regions may concentrate tiny grains or very 
fluffy aggregates (FA) (Wood 1998). TC may thus have been ubiquitous, helping to initiate 
the formation of "cometesimals" from porous grain aggregates at 10-30 AU or speeding the 
accretion of dust particles at high elevations (low gas densities) even in the terrestrial planet 
region. 

However, textural evidence might be difficult to obtain from this regime; subsequent 
compaction would obliterate evidence for any preferred size or density of easily squashed 
fluffy constituents. Chondrules and their parent chondrites, by nature of their availability 
and unique, persistent textures, provide the most obvious initial testing ground for turbulent 
concentration. 
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Figure 1: The optimal particle radius- density product rp s for turbulent concentration as a 
function of a at various locations in a typical solar nebula (converted from the optimal Stokes 
number St v = 1 using nominal nebula gas density and local rotation rate), (a) Minimum 
mass nebula {J- — 1); Solid silicate particles with chondrule sizes and densities (Ch) have the 
range rp s w 0.015 — 0.15, which would indicate the shaded range of a at 2.5 AU. Smaller gas 
densities at larger distances from the sun concentrate smaller rp s products (porous aggregates 
or PA) and even "fluffier" d=2 fractal aggregates (FA) or their monomers, for any a. (b) 
Somewhat different results obtained assuming an enhanced surface mass density JF = 10 
times higher than the minimum mass requirement. 



4 Concentrated particle size distribution 

If dense clusters of particles are precursors of primitive bodies in any way, the relative abun- 
dance as a function of Stokes number St v for particles in dense regions should show some 
correspondence to the distributions found in chondrites. To further compare the predictions 
of TC with meteorite evidence, we have recently determined the detailed form of the size 
distribution of selectively concentrated particles (Hogan and Cuzzi 2000). Numerical simu- 
lations were performed with homogeneous, isotropic, incompressible 3D turbulence at three 
Reynolds numbers Re = 62, 246, and 765 (these values differ from those given in Hogan et al. 
1999 because we now adopt the more correct definition of Vj, as one component of the turbu- 
lent velocity rather than its magnitude). The particles can be given arbitrary aerodynamic 
stopping times t s ; their motions respond only to gas drag and are integrated in the spatial 
domain. Feedback by the particles onto the gas is not incorporated. The computationally in- 
tensive calculations are run on 16 Cray C90 cpus at Ames Research Center. Simulations were 
initiated with uniform spatial distributions of particles, themselves uniformly distributed in 
stopping time over the range St v = 0.1-6, and the relative equilibrium abundance of par- 
ticles was studied as a function of St v and C. In the large-C limit of interest, the shape of 
the relative abundance distribution A(St rj ) was found to be independent of both C and Re 
and only very weakly dependent on the spatial binning assumed (relative to rj) (Hogan and 
Cuzzi 2000). Thus, we believe the numerical results should be valid as a prediction of the 
size distribution in dense particle concentrations under nebula conditions. The theoretical 
results are adequately fit by a lognormal distribution over the core range St v = 0.5 to 2.0. 

We compared our predictions with binned relative abundance data for chondrules disag- 
gregated from one carbonaceous and four ordinary chondrites (Paque and Cuzzi 1997, 2000; 
Cuzzi et al 1999). While a few ostensibly similar data sets exist in the literature (eg. Hughes 
1978, Eisenhour 1996, Rubin and Grossman 1987), they generally rely on radii measured 
microscopically from thin sections of meteorites, or even from disaggregated chondrules, and 
adopt some average chondrule density pj. However, recall from equation (6) that the aero- 
dynamic stopping time t s , which selects particles for TC or any other aerodynamic sorting 
process, depends on the product rp s for each object. Data for these meteorites imply that 
merely measuring the chondrule radius distribution and assuming some mean density will 
slightly, but noticeably, misrepresent the rp s distribution because of chondrule-to-chondrule 
density variations (Cuzzi et al 1999, Paque and Cuzzi 2000). 

In figure 2 we compare the particle abundance shape function A(St) from our numerical 
simulations (Hogan and Cuzzi 2000) with the meteorite data. The meteorite histograms were 
fit by lognormal functions and aligned horizontally with the predicted histogram by assuming 
that the optimally concentrated particle size-density product corresponds to St v = 1. We 
conclude that TC by itself can explain the very narrow chondrule size distribution, whatever 
the chondrule formation process may have produced. There is evidence that both larger and 
smaller chondrules were created; TC would predict that these non-optimally sized particles 
were simply not concentrated to a sufficient degree for eventual incorporation into meteorites 
unless they were swept up into chondrule rims (if originally smaller) or broken to the proper 
size (if originally larger). Some particles of arbitrary size will always appear just by accident, 
perhaps captured by a dense cluster; parent body processes will confuse the situation further. 
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Figure 2: Comparison of the size-density distributions of chondrules disaggregated from four 
ordinary chondrites (solid symbols; Paque and Cuzzi 2000) with the theoretically determined, 
.Re-independent, shape distribution for the relative abundance of particles concentrated in 
turbulence as a function of St v (open symbols; Hogan and Cuzzi 2000). 



10 



5 Probability distributions and multifractals: 

As mentioned earlier, our numerical results (CDH96) are obtained at Re which are six orders 
of magnitude smaller than plausible nebula Re. Thus, it is very important to understand the 
.Re-scaling of TC properties in detail. We believe this has become possible using a connection 
between TC and the .Re-independent scaling properties of fractals, which are rooted in the 
properties of the turbulent cascade itself (Hogan et al 1999; henceforth HCD99). Fractals and 
multifractals are readily associated with cascade processes. Turbulence (and specifically its 
inertial range) is the archetype of a cascade process (Tennekes and Lumley 1972, Meneveau 
and Sreenivasan 1987, 1991), and considerable attention has been devoted to multifractals 
in turbulence (Frisch and Parisi 1985, Mandelbrot 1989). 

A working definition of a fractal is a structure which is generated by sequential application 
of a scale-invariant rule on regularly decreasing spatial scales. Simple fractals with constant 
(but non-integer) dimension are invariant to changes in scale, and result from rules which 
produce a binary distribution of the local density (say, either or 1). Examples of these 
are the Cantor set or the Sierpinski gasket (Mandelbrot 1989), in which segments of a line, 
or portions of a plane, are simply removed without changing the surrounding values. Their 
average density, as a function of scale e, may be written as p(e) = po e ~ d with dimension 
d = 0.63 (between a point and a line) and d = 1.52 (between a line and a plane) respectively. 

In contrast, multifractals result from application of rules (or probability distributions of 
rules) in which the local measure is changed, while conserving the total measure, by unequal 
repartitioning of the content of a bin into sub-bins of regularly decreasing spatial scale. For 
example, it is easily seen that partitioning some quantity into the two halves of a bin with 
unequal proportions (say 0.7 and 0.3) raises the mean density in the first half and decreases 
it in the second. Repetition of this rule produces some bins which become ever denser with 
decreasing scale, and others which become ever less dense - with all combinations in between. 
The ensuing spatial distribution has no well defined local value in the limit of diminishing 
bin size; that is, the local values are spatially spiky, "intermittent", or "singular" (Mene- 
veau and Sreenivasan 1987 present a short and readable discussion of cascade processees 
and multifractals). The spatial distributions of multifractals are, however, predictable in a 
statistical sense, using Probability Distribution Functions (PDFs) which are derived directly 
from their dimensions. For instance, dissipation of turbulent kinetic energy, which occurs on 
the Kolmogorov spatial scale, is not spatially uniform but has the spatial distribution of a 
mutifractal (Chhabra et al 1989). 

In multifractals, the dimension varies with the value of the measure. In our case, the 
measure is particle concentration factor C, which is the ratio of the local particle volume 
density to its global average value. The fractional probability Pi of a particle lying in a bin 
which contains iVj particles out of N p total particles is defined as Pi = Ni/N p = e ai , where 
e = bin size/domain size. The scaling index a can be viewed as a local dimension for P. 
The associated concentration factor Cj = J, , where u- is the volume of a bin and v is the 

total volume of the domain; thus Ci = *', p = (e ai /e 3 ) = e ai ~ 3 . Expressing the bin size, or 
scale, as some multiple B of r\, and the domain size in units of the integral scale as DL, we 
find the domain normalized bin size to be e = Br]/ DL = -j|.Re _3//4 = 1/TZ, where we have 
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Figure 3: Singularity spectra f(a) for particle concentration (solid symbols), for bin sizes 
of 2, 4, and 8 i], along with error bars representing deviations across many snapshots. The 
similar looking singularity spectrum for dissipation (open symbols) is from Meneveau and 
Sreenivasan (1987). Also shown (dashed curve) is our annealed f(a) for 2r] binning, covering 
C values which occur less than once per computational domain (and thus extends to negative 
values). Figure adapted from HCD99. 



also used the inertial range relationship r/ = LRe 3 / 4 . Thus, C = TZ 3 a , or 

InC 

a = o — — — -. 
InTZ 



(9) 



The normalized PDF for a is usually written as a fractional volume F v (a) = p(a)e 3 ~^ a \ where 
the important function f(a), called the singularity spectrum (Halsey et al 1986, Chhabra et 
al 1989, Mandelbrot 1989), plays the role of a dimension for F(a). The prefactor function 
p(a) is only weakly dependent on scale, and can be approximated as y/lnR, (Chhabra et al 
1989). The function f(a) is discussed more in the next subsection. 

We define the PDF F V (C) as the volume fraction occupied by bins having concentration 
factor within (C, C + dC), with l^ZZ F v(C)dC = J^ F(a)da = 1. Transforming variables 
and their PDFs, we get 



F V (C) = F{a) 



da 

dC 



p(a)K f ^- 3 TZ f ^~ 



(10) 



Clnft CVhtfl' 

The fraction of particles occupying bins within (C, C+dC) is defined as F p (C) = CF V (C), 
and both F V {C) and F P (C) have cumulative versions F v (> C) = J^ 1 F v (C)dC, and F p (> 
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C) = f£?F p (C)dC. For a stationary, ergodic process, one expects that the fraction of 
particles having C at any given time (F p (C)) is identical to the fraction of time spent in 
regions of concentration C by a typical particle, defined as F t (C). Below, this is confirmed 
numerically. 

5.1 Numerical results at low Re: 

5.1.1 PDFs for concentration 

Using the methodology of Chhabra et al (1989), we have shown that the spatial distribution 
of optimally concentrated (St v = 1) particles is a multifractal, and that its singularity 
spectrum f(a) is invariant over more than an order of magnitude from Re = 62 to Re = 765 
(HCD99, who actually give Taylor microscale Reynolds numbers Re\ = \/l5i?e = 40 to 
140). This i?e-invariance for f(a) is seen only when binning is done on some fundamental 
flow-relative scale such as a multiple of the Kolmogorov scale rj. The function f(a) is found 
by applying weighted binning techniques to the distribution of interest (Chhabra et al 1989, 
HCD99), and depends on both D and B - that is, on the number of integral scales sampled 
as well as on the binning scale. Careful reassessment of different definitions of the integral 
scale (see Hinze 1975, chapter 2) has led us to values slightly smaller than those published by 
HCD99, implying that the computational domain appropriate for our f(a) is 9L on a side. 
Our results for f(a) relative to this domain are shown in Figure 3 (HCD99). The sets with 
error bars are "quenched", or averaged over samples each DL in extent. The fact that there 
is no less than one cell in each sample found at maximum concentration C max = 7£ 3 ~ amin 
results in /(a m i n ) = for these sets. This can be seen by setting the product of either F v (a) 
or F V (C) by the normalized domain volume 1Z Z equal to unity, after allowing for the fact 
that the F v are differential functions of their arguments by multiplying them by a m { n or C max 
respectively. 

Information about higher (less probable) concentrations than are typically found in a 
volume DL on a side must come from analyzing a large number of realizations of each sample 
and determining an "annealed" average, where /(a m i n ) < (eg., Chhabra and Sreenivasan 
1991). The annealed version (dashed curve in figure 3), also binned by 2r], is well fit by 
the function f(a) = -12.414 + 89.659/a - 132.01/a 2 , for a > a min = 2.0. We have restricted 
ourselves to 2i] binning for the present to preserve good statistics; smaller binning scales 
sample a deeper cascade and will generate large C values with higher probability. The three 
sets of f(a) in Figure 3, as binned over 2, 4, and 8 rj illustrate the multifractal or "singular" 
nature of the distribution. For any given Re, smaller a values correspond to larger C values 
(equation 9), and larger values of f(a) imply larger F v . The smaller a m i n , or larger f(a) 
for a > a m i n seen for the smaller binning scales implies that, averaged over a bin, smaller 
binning scales retain large C far more commonly than larger binning scales. That is, as the 
scale is reduced, no asymptotic or well-defined limiting local value is reached. Also shown 
in Figure 3 is the singularity spectrum for dissipation of turbulent kinetic energy (open 
symbols), which is scale-independent. The agreement between this spectrum and that for 
particle density binned at the 2i] scale is intriguing. The deviation between our spectra 
(solid symbols and dashed line) and that for dissipation seen toward the right hand side 
(large a) is due to incomplete sampling of very low particle density regions in our particle 
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Concentration factor C 

Figure 4: Probability Distribution Functions (PDF's) for the fraction of particles lying in 
regions with concentration factor greater than C {F p (> C)), or, equivalently, the fraction of 
time spent by any particle in such regions (F t (> C)). The three sets of points are binned 
directly from our numerical simulations; the associated curves are calculated from a single 
averaged f(a) obtained from all three values of Re (dashed curve in Figure 3). As discussed 
in section 5.2, the curves without points use the same f(a) to predict PDF's at the larger Re 
corresponding to four plausible nebula a values: 10~ 2 (solid line), 10~ 3 (long dashed line), 
10~ 4 (short dashed line), and 10~ 5 (dotted line). 

density simulations because of memory limitations on the number of particles we can follow. 
Dissipation, being a continuously varying function, is not subject to this effect. For this 
reason, F v (> C) was less well defined than F p (> C), which, by definition as a particle rather 
than volume fraction is always fully characterized. In any case, our primary interest is in 
the zones of high concentration (small a). 

As a check on the method (subtle normalization issues are discussed in Chhabra et al 
1989), an average f(a), obtained from our calculations at all three Re values was used 
to calculate F p (> C) at each of our three Re values, for comparison with the F p (> C) 
distributions determined directly from numerical results (HCD99). As shown in figure 4, 
the single f(a) does quite well at predicting F p (> C) at all three Re, even though they 
all subtended slightly different numbers of integral scales. For any C, F p (> C) increases 
with Re and thus R. This may be puzzling after only a quick inspection of equation (10), 
since f(a) — 3 is always negative (figure 3). However, for a fixed C, a(C) increases with 
1Z (equation 9), and thus f(a) also increases (figure 3), so the exponent in equation (10) 
becomes less negative. Since any function is more sensitive to its exponent (here f(a) — 3) 
than to its base (here 1Z), equation (10) implies F v (> C) and thus F p (> C) increase with 1Z. 
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5.1.2 Particle time histories 

The particle time histories are of interest in their own right, and illustrate how particles 
experience a fluctuating background concentration as they wander through the fluid. Figure 
5 illustrates histories for several randomly chosen particles. The particles are moving at 
roughly constant space velocity (approximately the velocity Vl of the largest eddies, since 
they are trapped to nearly all eddies), and repeatedly encounter zones of different density 
with little noticeable effect on their velocity (their stopping time is much longer than the 
clump transit time tj/Vl)- In these simulations, where no particle interactions are computed, 
the particles pass through the dense zones and continue their evolution. The dense zones per 
se persist for times much longer than the passage time of a single particle (CDH96, She et 
al 1990). The more numerous, lower density zones are encountered more frequently, and the 
rare, very high density zones less frequently. This time history, essentially a (convoluted) 
ID path through the computational volume, has the same "intermittent" or spiky structure 
as seen for dissipation (see, eg., Meneveau and Sreenivasan 1987; Chhabra et al 1989). 

From simulations such as these, we have validated the ergodic assumption that the frac- 
tion F p (> C) (spatially averaged over all particles at several snapshots in time) is equal to 
the fraction of time F t (> C) spent by a given particle in regions denser than C (temporally 
averaged over extended trajectories for a few particles). The comparison is shown in figure 
6. A "random walk" calculation does not accurately reflect this behavior; the particles are 
not moving randomly through the volume and encountering dense zones with probability 
given by their volume fraction; rather, particle trajectories are "attracted" to the dense 
zones, or, probably more physically, repelled from the sparse, complementary eddy zones. 
This behavior is correctly captured using F P (C). 

5.2 Predictions for nebula Re: 

The properties of multifractals make it much easier for us to predict nebula properties than 
merely making extrapolations from Re values we have studied numerically. In the multifrac- 
tal context, equation (10), and the Re- independence of /(a), suggest identification of f(a) 
with the level-independent "rule" , and 1Z oc Ljr\ = Re 3 ' 4 " with the number of steps in the 
cascade. As Re increases, the inertial range (or number of eddy bifurcations) between L and 
7] also increases. The cascade model described by Meneveau and Sreenivasan (1987), men- 
tioned at the beginning of this section, is a specific case of the binomial cascade discussed by 
Mandelbrot (1989); all provide insight as to how higher C result from a deeper cascade with 
more steps. Mathematically, equations (9-10) show how high C values become more likely. 
If f(a) is indeed a i?e-independent, universal function for optimally concentrated (St v = 1) 
particles, we can then predict PDF's for any Re (i.e., any given nebula a) using equation 
(10). 

In addition to our own experiments over an order of magnitude in Re, several other lines 
of argument support the assumption of a level-independent and .Re-independent "rule" . Dis- 
sipation of k and particle concentration are physically connected through their mutual pref- 
erence for the Kolmogorov scale; also, the shapes of the singularity spectra for dissipation, 
and for particle concentration binned at close to the Kolmogorov scale, are similar (figure 3; 
HCD99). The singularity spectrum of dissipation has already been connected to the turbulent 
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Figure 5: Ambient concentration encountered by four different particles moving in 3D turbu- 
lence, as a function of time (measured in large eddy turnover times LET = uj(L)~ v ), without 
feedback onto the gas. These may also be regarded as a longer history for a single particle. 
Note how denser regions are encountered less frequently. 
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Figure 6: Comparison of F p (> C) from binned data in several snapshots (squares, temporal 
realizations) containing all particles (as in Figure 4), with F t (> C) (triangles), calculated 
from extended time histories of 16 different particles (as in Figure 5). 
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cascade process (Meneveau and Sreenivasan 1987). Thus, turbulent concentration is proba- 
bly also closely related to the turbulent cascade process. The turbulent cascade is known to 
have Re- independent properties in the inertial range. For instance, the f(a) for dissipation 
has been shown to be .Re-independent - from numerical work at Re ~ 100, including our own, 
through laboratory experiments with Re ~ 10 4 , to experimental studies of the atmospheric 
boundary layer with Re ~ 10 6 (Chhabra et al 1989, Hosogawa and Yamamoto 1990). Fur- 
thermore, it has been directly demonstrated from analysis of observations of dissipation in 
large Re turbulence that the probability distribution of "partition factors" or "multipliers" 
is independent of level within the inertial range (Sreenivasan and Stolovitsky 1995). Based 
on these arguments, we believe and presume that the particle concentration singularity spec- 
trum f(a) remains Re- independent in turbulent cascades with far larger values of Re than 
those of our numerical experiments. Given this invariance, we can predict nebula conditions 
more confidently than from extrapolation alone, as was done by CDH96. We obtain PDF's 
at Re which are much larger than directly accessible values by containing the .Re-dependence 
purely within 1Z = DL/ Br] = (D/B)Re 3 ^. As argued at the beginning of this section, the 
f(a) is associated with a certain D and B, which we interpret as being the D m 9 and B=2 
implicit in our f(a) (section 5.1.1). The nebula TZ is then only dependent on Re and easily 
determined for any turbulent a as 1Z = 1/e = DL/ Br] = (D/B)Re 3 ^ 4 ~ 4.5i?e 3 / 4 . Recalling 
v m = mH 2 c/2p g <TH 2 (section 3), 

* = 4.5 f^f = 4.5 MM 3 " ^7 x 10= f^J— f ( °)«' . (11) 

Also shown in figure 4 are predictions of F p (> C) = F t (> C) for four values of nebula a, 
using equations (10) and (11). Based on these predictions, St v = 1 particles spend about 
1 - 10% of their time in regions with C > 10 3 under nebula conditions at different a. All 
curves in figure 4 assume a minimum mass nebula (£ = 430 g cm 2 at 2.5 AU), but mass 
enhancement by some factor F would play a role similar to a (see also section 3). 

6 Implications and discussion 

6.1 Regions of moderately high density 

For concentration factors leading to p c h/Pg at least as large as unity, particle mass loading 
might start to affect the "partition rules" and change the statistics of the cascade. Neverthe- 
less, several interesting effects can result from concentrations less than this possible limit. If 
the region in question has an average mass density in solids of p p = f so \Pg = 5 — 8 x 10~ 3 p 9 
(at, say, 400°K; Pollack et al 1994), and some fraction / c h < 1 of this amount is in chondrule- 
sized particles (recalling that other particle sizes are not susceptible to TC), it is clear that 
C as large as (fsoifch)" 1 , at least several hundred (for / c h = 1), remains free of mass loading 
concerns. If / c h << 1, far larger C are allowed. The point is that, depending on / so i/ c h, C 
can get quite large without violating the mass loading caveat. 

While most of the volume of the nebula is characterized by C < 100, Figure 4 shows 
that chondrules in the nebula spend a significant fraction of their time residing in regions 
where the particle density is enhanced by a factor of C > 100. These results might help us 



understand why some chondrule types seem to have been heated under unusually oxidizing 
conditions (eg., Rubin et al 1989), which has been attributed to vaporization of a large 
ambient density of solids - enhanced by several orders of magnitude over solar - as part of 
the chondrule formation process (note, however, that it is not universally accepted that the 
chondrule oxidation state information necessarily implies a large background abundance of 
oxygen; Grossman 1989). 

Several chondrule melting processes - shock waves (Hood and Kring 1996, Connolly and 
Love 1998) or lightning bolts (Desch and Cuzzi 2000, Desch 2000) are envisioned to occur 
ubiquitously within the region of the nebula in which we propose TC operates to sort the 
chondrules after their formation. In fact, Desch and Cuzzi (2000) showed that TC itself is 
an enabling factor in generation of nebula lightning. They found that the optimal conditions 
for energetic lightning bolts are found in 1000-km regions with C ~ 100 - especially if the 
nebula were denser than "minimum mass" - T ~ 10. Either lightning or shock waves - 
heating dense zones where the background density of solids is large (cf. Hood and Horanyi 
1993, Hood and Kring 1996, Connolly and Love 1998) - could elevate the local Oxygen 
abundance by evaporating some part of the solids. However, TC does not concentrate "fine 
dust" - only chondrule-sized particles - so significant local volatilization of some chondrules, 
or the surfaces thereof, would be implied. Whether this is consistent with mineralogical 
signatures which are seen in the survivors would be useful to address in the future. 

6.2 Encounters with very dense regions 

The scenario of CDH96 suggested in a qualitative way that TC could lead to stable, if low- 
density, clusters or clumps of chondrules as direct precursors of planetesimals. We delineate 
the logic of this argument, using the PDF's derived in section 5, and then discuss the 
difficulties with the original scenario. 

The PDF's can be combined with the particle velocity through space V p to calculate 
the "encounter time" T cnc of a chondrule with a region of arbitrary C, using a duty cycle 
argument. The fraction of time spent by a typical chondrule-size particle in clumps with 
concentration greater than C, F t (> C) « ^p ~ ^T ' wnere M> c ) « T cnc is the 
time spent traversing a bin with concentration greater than C (section 5.1.2). Thus, for bins 
of dimension 2r\ and particle velocity V p , t m (> C) = 2r//V p , so T cnc ps ^"[^J = f(>c) = 
y- F } >c) , and the encounter rate is jr— = -f-F p (> C). We have verified numerically that, 
as expected for particles with stopping times t s much shorter than the overturn time of the 
largest eddies f^" 1 , V p is nearly identical to the typical turbulent gas velocity \/2k ?s c^fci 
(Volk et al 1980, Markievicz et al 1991). 

We may calculate the encounter rate (and time) with a cloud so dense that a particle 
becomes entrapped with its neighbors and possibly removed from further free circulation. 
Normally, as seen in figure 5, particles traverse dense regions without incident, because their 
gas drag stopping time is longer than their transit time (CDH96). An entrapment threshold 
occurs if interparticle collisions can prevent particles from passing through a cloud; this 
implies a critical cloud optical depth r coll of unity, defining a critical C co \\. For small, dense 
2??-sized clumps, r co n = 1 = 2r]C con (p ch /m)7i(2r) 2 , where p ch = f so \f ch p g is the average 
(unconcentrated) chondrule mass density, / so i ~ 5 x 10~ 3 at 400K is the fractional mass in 
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solids, / c h is the fraction of solids in chondrules, and m and r are chondrule mass and radius. 

Then 

4x 10 5 /lO-^X 174 
C co n w 33rp s /r) Pg f ch « — . (12) 

/ch V " / 

In the final expression above, we have substituted the first expression of equation (8) for the 
product rp s , since the value of this product for optimally concentrated particles is constrained 
by nebula conditions (section 3). We substituted the general expression for r\ = LRe 3 ^ 4 = 
(H J a) l / 4 {v m / c) 3 / 4 (~ 0.5km (lO^/a) 1 / 4 ), and assumed other nominal parameters at 2.5 AU 
(£ = 430 g cm -2 , c = 1.5 x 10 5 cm/sec). 

For nominal, minimum mass nebula parameters, assuming / c h ~ 1, and using the PDF's 
we have in hand (binned at 27/ scales), we find that T cnc varies between 10 3 years for a = 10~ 2 
to 10 8 years for a = 10 -4 . The encounter times are sensitive to nebula parameters adopted 
and would decrease for smaller binning (eg., by B — 1), or enhanced nebula densities [T > 1). 
This quantifies the length of time chondrules spend freely wandering before encountering a 
region in which they are certain to undergo collisions. We suggested earlier that such an 
encounter at T enc (C co n) removes the chondrule from further circulation by entrapping it with 
others in the dense cluster. Alternately, such dense zones might merely provide a large in- 
crease in the collisional aggregation rate of optimally sized particles. In fact, if chondrules 
or groups of chondrules undergo collisional aggregation at low relative velocities (and par- 
ticles of comparable size always have low relative velocities in turbulence; Weidenschilling 
and Cuzzi 1993), they may form fractal aggregates with mass proportional to radius squared 
(Weidenschilling and Cuzzi 1993, Beckwith et al 2000). The average density of such aggre- 
gates is inversely proportional to their bounding "radius" , so they retain the same stopping 
time as their individual components (individual chondrules) , and can continue to participate 
inTC. 

However, before this line of thought can be pursued much further, mass loading (discussed 
in the next section) must be assessed. Depending on their linear extent, the optically thick 
clouds described above can reach a mass density orders of magnitude larger than that of the 
gas, which would probably invalidate our assumption of no feedback. In general, for a clump 
of mass density p c h composed of chondrules of radius r and density p s , 

rpsT C o\\ pch rp s 

Pch = QD ; or — = — for Ten = 1. (13) 

Wr] p g 2>B-qp g 

In deriving equation (13) we have expressed the linear size of a cluster as Brj. For rj ~ 0.5 
km (lO^/a) 1 / 4 and p g = 1.1 x 10~ 10 g cm -3 , the particle density exceeds 10p g for r co ii = 1 
clusters smaller than (60 — 600/jF) km - about (100 — 1000/jF)r/ for rp s in the normal range 
for chondrules ~ 0.02 — 0.02 (figure 1). 

6.3 Limitations due to high mass loading: 

The cascade process model is only valid as long as no new physics emerges at some step in 
the cascade to change the partition factors, as represented globally by f(a), applicable to 
subsequent steps in the cascade. However, application of the expressions above for optically 
thick clusters of scale 2rj implies p c h/p g ~ 3000. While these high concentrations literally 
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relate to regions comparable to 77 in size, where there is no turbulence to damp, the cascade 
that is needed to produce such a cluster must have extended to larger sizes in the surrounding 
"penumbra" where the particle density, while lower, might still be large enough to damp tur- 
bulent motions. We have made some preliminary calculations of f(a) from low- .Re numerical 
simulations having mass loading p p /p g ~ 1. Even here, turbulent concentration persists (to 
p p ~ 30p 9 ), but the f(a) is altered in the sense that high C values have a lower probability. 
For comparison, unloaded simulations at this Re result in C max ~ 60. The magnitude of the 
mass loading effect is approximately equal and opposite to that of decreasing the bin size 
from 2i] to 77. Clearly, this effect must be better quantified before more specific predictions 
of T enc and chondrule accumulation timescales can be made. Even for small a ~ 10~ 3 , the 
mass loading regime of questionable validity is not a large one (note that equation (13) above 
showed that regions of size 30077, or 100-300 km, are probably within the range of validity); 
however, it is certainly an interesting one. If we further restrict our attention to optimally 
concentrated particles by combining equations (7 or 8) and (13), we obtain 

p ±= (j^_Y 2 (^_\(h_y 2 (14) 

l>, V4V2ff H J \3V2Bii) \ap,) 

which shows that regions which are denser and/or have more intense turbulence (higher a), 
such as may have characterized the very early evolutionary stages of the nebula, are less 
prone to mass loading difficulties. That is, they provide r co ii = 1 at lower p c h/p g - Figure 
7 shows that conditions at 1AU (gas density p g ~ 10~ 9 g cm -3 ), with a ~ 10~ 2 — 10 _1 , 
provide r co n = 1 with p p /p g < 100. Full treatment of a wider parameter space, allowance 
for mass loading and other binning scales B, and connection to hypotheses for chondrule 
formation (Desch and Cuzzi 2000) and dust rimming (Cuzzi et al 1999, Morfill and Durisen 
1998) which are tied to TC, will be treated in subsequent publications. 

7 Discussion: 

Turbulent concentration promises to provide another tool for our planetesimal construction 
toolkit. Nominal nebula conditions, including weak turbulence, lead to concentration fac- 
tors large enough to concentrate chondrule- sized particles to densities greatly exceeding the 
nominal nebula average. The presence of such dense regions might help us understand some 
aspects of chondrule formation. Generalization of the physics shows that fluffy aggregates 
would also be subject to turbulent concentration, but in nebula regions where the gas den- 
sity is lower and/or the turbulent intensity higher than expected for the terrestrial planet 
region. If TC is a key aspect of the primary accretion of planetesimals, any chondrules 
formed or transported to low density, outer solar system regimes might only have escaped 
incorporation into planetesimals to a greater degree than in the terrestrial planet region. 
The size distributions produced by turbulent concentration are in very good agreement with 
those found in primitive chondrites; however, the fact that most primitive meteorites con- 
tain ample evidence for abrasion, fragmentation, and other mechanical processes which may 
well have continued long after the nebula gas vanished and aerodynamic processes became 
irrelevant, warns us that we should not expect any single process to explain all properties of 
even "primitive" meteorites. 
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Figure 7: Variation of p p /p g with a (for cluster optical depth unity) for the same nebula 
regions as shown in figure 1. Regions with more severe turbulence (larger a) select smaller 
particles for concentration which, having larger area per unit mass, can be concentrated to 
T"coii without reaching such large mass loading ratios as under more "nominal" conditions as 
shown in figure 1. 
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The spatial distributions of concentrated chondrules under nebula conditions can be 
predicted using the multifractal properties of turbulent concentration. Various timescales 
of interest can then be estimated, including the timescale in which a chondrule encounters 
a cluster so dense that it is entrapped. However, preliminary estimates of timescales for 
simple accumulation by entrapment in unusually dense clusters are comfortably shorter than 
nebula lifetimes only for higher a values than those which best explain chondrule size-density 
products. 

Mass loading is expected to change the nature of the process at some point in the cascade 
where the concentrated particle density becomes comparable to, or greater than, that of the 
gas. Depending on the mass fraction in solid "chondrules" , this might imply concentration 
factors as small as 10 2 - or possibly far larger. Mass loading has not yet been quantitatively 
treated, but it seems it will inevitably increase the encounter timescales by decreasing the 
abundance of very dense clusters. Perhaps turbulent concentration only augments collisional 
accumulation for optimally sized particles. 

The evolution of dense clusters in the presence of mass loading and the vertical component 
of solar gravity deserves study (eg. Wang and Maxey 1993). Do clusters retain their identity 
or disperse? Would they be compacted or dispersed by settling, by collisions with other dense 
clusters en route to, or within the midplane, or by shock waves in the nebula? Perhaps only 
clusters formed near the midplane avoid the dispersion due to settling. 
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